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Multi-boson symmetrization effects on two-particle Bose- 
Einstein interferometry are studied for ensembles with arbi- 
trary multiplicity distributions. This generalizes the previ- 
ously studied case of a Poissonian input multiplicity distri- 
bution. In the general case we find interesting residual cor- 
relations which require a modified framework for extracting 
information on the source geometry from two-particle corre- 
lation measurements. In sources with high phase-space den- 
sities, multi-boson effects modify the Hanbury Brown-Twiss 
(HBT) radius parameters and simultaneously generate strong 
residual correlations. We clarify their effect on the correla- 
tion strength (intercept parameter) and thus explain a variety 
of previously reported puzzling multi-boson symmetrization 
phenomena. Using a class of analytically solvable Gaussian 
source models, with and without space-momentum correla- 
tions, we present a comprehensive overview of multi-boson 
symmetrization effects on particle interferometry. For event 
ensembles of (approximately) fixed multiplicity, the residual 
correlations lead to a minimum in the correlation function 
at non-zero relative momentum, which can be practically ex- 
ploited to search, in a model-independent way, for multi-boson 
symmetrization effects in high-energy heavy-ion experiments. 
PACS number(s): 25.70-q, 25.70Pq, 25.70Gh 



I. INTRODUCTION 

In high energy hadronic and nuclear collisions large 
numbers of pions are created (up to several thousand in 
Au+Au or Pb+Pb collisions at the Relativistic Heavy Ion 
Collider (RHIC) or the Large Hadron Collider (LHC)). 
This raises the issue of multi-particle symmetrization ef- 
fects among the final state particles and their effects on 
the particle spectra and momentum-space correlations. 
If the particles are set free from a state of sufficiently 
high phase-space density, such effects are expected to be- 
come measurable. This has triggered many investigations 
of multi-boson symmetrization effects over the years (see 
p|-p^| and further references given in those papers). In 
particular the possibility of Bose-Einstein condensation 
of pions in heavy-ion collisions (sometimes misleading! 



called "pion laser" ) has attracted attention 
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Recent analyses of experimental data point, however, to 
universally low pion phase-space densities at freeze-out 
|§HT§. On the other hand, the authors of Ref. (Dp 
present a chemical analysis of the hadronic final state in 
158 A GeV/c Pb+Pb collisions at the CERN Super Pro- 
ton Synchrotron (SPS) which yields a large pion chemical 
potential, fj, n ~ m^, corresponding to rather large pion 
phase-space densities. The theoretical investigation of 
multi-boson symmetrization effects thus continues to be 
a challenging problem of practical relevance. 

In this paper we present a comprehensive analysis of 
multi-boson symmetrization effects on multiparticle spec- 
tra and Bose-Einstein interferometry. Our work not only 
provides a synthesis of many previously published investi- 
gations in a common and very general formal framework, 
but also extends them in one important aspect: while 
nearly all previous studies (the only exception known to 
us is Ref. [plj|) analyzed either event ensembles of fixed 
particle multiplicity or an ensemble, in which multi-boson 
symmetrization effects at finite phase -space density mod- 
ified in a very specific way (see Sec. [II C) a Poissonian 
input multiplicity distribution at vanishing phase-space 
density our formalism allows to take into account 
arbitrary measured multiplicity distributions. In addi- 
tion to the now well-established "Bose clustering effect" , 
by which the effective emission function becomes nar- 
rower in both coordinate and momentum space, we show 
that multiboson symmetrization in general also intro- 
duces "residual correlations" in the two-particle corre- 
lation function, i.e. deviations from unity which, in con- 
trast to the Hanbury Brown-Twiss (HBT) effect result- 
ing from two-particle exchange, are not directly related 
to the source size. These residual correlations modify the 
effective correlation strength A, i.e. the ratio of the two- 
particle correlation function at zero and infinite relative 
momentum, and they introduce a new dependence of the 
correlator on the relative momentum q which modifies 
the extraction of the source size via HBT radii. 

These residual correlations are generic. We know only 
two multiplicity distributions for which they do not arise, 
the "Poisson limit" discussed in Sec. [IIIC| 
and the "covariant current formalism" discussed in Ap- 
pendix |A] {200. In both cases, one prescribes an "in- 
put" multiplicity distribution (which refers to the infinite 
phase-space volume limit where multi-boson symmetriza- 
tion effects can be neglected Jl9)), whereas the actually 
measured multiplicity distribution is a complicated func- 
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tion of the average phase-space density of the source. On 
the other hand, residual correlations do arise, in particu- 
lar, for events with fixed multiplicity; there they lead to 
strong effects in sources with large phase-space density 
which were noted before [jl]j6||| but whose origin has up 
to now not been fully clarified. 

The residual correlations and Bose clustering effects 
disappear for dilute sources which emit particles from 
a state of low phase-space density. Only in this limit 
the usual formalism (see, for example, the recent reviews 
p2| ) for extracting effective source sizes from two-particle 
correlation measurements can be applied. As the source 
phase-space density (more accurately: the spatially aver- 
aged phase-space density (/) (K) |l3|-|l(| of particles with 
pair momentum K at which the correlation function is 
evaluated) increases beyond a certain threshold, multi- 
boson effects set in rather abruptly. When this happens, 
a more general formalism must be applied to separate 
residual from HBT correlations. This is worked out in 
the paper. 

An important stepping stone for the discussion of 
multi-boson symmetrization effects was provided 2 years 
ago by Csorgo and Zimanyi H who found an analytic 
solution for the corresponding recursion relations within 
a class of Gaussian source models. This analytic solu- 
tion dramatically reduces the complexity of the problem 
and has thus triggered many recent (re-)investigations of 
multi-boson effects; our paper is no exception. Before us- 
ing them we rewrite their results in a physically more 
instructive way which clearly exhibits the relevant phys- 
ical parameters in the problem. This is straightforward, 
but we have not seen it published before, and it has im- 
portant practical consequences. In the applications we 
focus on Bose-Einstein interferometry studies of the two- 
particle correlation function; a detailed analysis of multi- 
boson effects on the measured single-particle spectra was 
recently published in fl^| . We argue that, once the effects 
of multi-boson symmetrization on the normalization and 
shape of the two-particle correlation function are qualita- 
tively understood, they can be efficiently used to search 
experimentally for multi-boson effects and for evidence 
of large phase-space densities in high energy heavy ion 
collisions. This may help to settle the controversy over 
whether or not pions are emitted with large chemical po- 
tentials 

Our paper is organized as follows: In Section || we 
write down a very general form of the density matrix 
for the ensemble of pion emitting sources in terms of a 
superposition of elementary source currents. The lat- 
ter are assumed to have random phases to ensure inde- 
pendent or "chaotic" particle emission. We introduce 
the source Wigner density and the "ring integrals" ||] 
from which the normalization of the density operator 
and later the particle spectra can be calculated. Care- 
fully normalized expressions for the latter are presented 
in Section [II , both for events with fixed multiplicity and 
for multiplicity-averaged events. We define the correla- 
tion strength and the incoherence parameter (the latter 



equals 1 for the fully chaotic sources studied in this pa- 
per), and we discuss specifically the previously studied 
"Poisson limit" in which for finite phase-space volumes 
multi-boson effects modify in a specific way a Poissonian 
input multiplicity distribution at infinite phase-space vol- 
ume. In Section IV] we review and reformulate the ana- 
lytical solution p] of the multi-boson recurrence relations 
for Gaussian sources and show in particular that it can 
also be applied to Gaussian sources with x-p-correlations 
(the Zajc model |23|]) for which the momenta of the emit- 
ted particles depend on the emission point in space-time. 

At this point the basic formalism is complete, and 
in Section ^ we proceed to discuss numerical results. 
We separate the two-particle correlation function into 
Bose-Einstein and residual correlations and study the 
asymptotic normalization of the correlator, the correla- 
tion strength, the HBT radii and the range of the residual 
correlations as a function of the pion multiplicity distri- 
bution, the pion phase-space density, and the pair mo- 
mentum, covering essentially the complete physical pa- 
rameter space. The numerical calculations are checked 
against analytic results in the limits of large (v — > oo) and 
small phase-space volume (y — > I). For v — 1 the source 
saturates the quantum mechanical uncertainty limit, i.e. 
all pions occupy the same quantum state, forming a Bose 
condensate. This interesting limit is under complete ana- 
lytical control. Experimental consequences are discussed 
in Section VI: we argue that multi-boson effects can be 
detected by triggering on fixed multiplicity and looking 
for a minimum at non-zero relative momentum q in the 
corr elation function. Some concluding remarks in Section 
VII close the paper. 

In three Appendices we provide technical details: Ap- 
pendix |A| discusses how our formalism is related to the 
"covariant current formalism" of Gyulassy et al. pO| , pT| ], 
the only other known ensemble which does not suffer from 
residual correlations in the presence of multiboson effects. 
Appendices ^ and |^ provide formulae for the analytic 
calculation of residual and Bose-Einstein correlations in 
certain limits. 



II. SOURCE PARAMETRIZATION 

In this section we describe the parametrization of the 
density operator on which our analysis of multi-boson ef- 
fects will be based. Very little of the material in this 
section is original; the most important generalization of 
previous work, which almost exclusively discussed cither 
fixed pion multiplicity or the case of a Poissonian mul- 
tiplicity distribution in the absence of multi-bo son sy m- 
metrization effects (the "Poisson limit" of Sec. [II C ), is 
that we allow for an arbitrary pion multiplicity distribu- 
tion. Still, although many of the relations given in this 
section can be found scattered in the literature, we find 
it useful to summarize them and thereby establish our 
notation and our normalization conventions. This last 
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point is important since much confusion over multipar- 
ticle symmetrization effects has arisen from insufficient 
attention to normalization issues. 



A. The density operator 

With the normalized (measured) pion multiplicity dis- 
tribution p n we can write the density operator of the 
ensemble of collision events as 



oo 
n=0 



X>„ = 1, trp n = l, (2.1) 



n=0 



where p n is a normalized density operator for a system of 
exactly n pions. We generalize the formalism developed 
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|2J] and consider the ensemble as a superpo- 
sition of localized source amplitudes jo whose number N 
is distributed according to a normalized probability dis- 
tribution Pn (in general not related to p n above). The 
sources have identical structure, but they are centered 
at different points Xi and move with different 4- velocities 
Ui = 7i(l,Uj). Their average positions and velocities are 
distributed by a normalized classical phase-space distri- 
bution p(xi,Vi): 



d 4 Xi 7^ cPvi p(xi,v.i) = 1 



(2.2) 



Here we introduced the shorthand notation Q (with the 
Lorentz invariant measure dQ) for the phase-space co- 
ordinates (xi,Vi) of the sources, (-ff d 3 Vi = d 3 Ui/u° is 
the invariant measure for the boost velocities.) Further- 
more the source amplitudes jo are assumed to have a 
randomly distributed phase fa. We can thus write for 
the total source I21IJ2I 



A" 



J{x) = Y,j* i 3o{ki-{x-Xi)) 



(2.3) 



Following d;@,|2C||l],|Il we make tne ansatz 



with 



|n) = |n[JV; {&&}]> = -= 



n)(n\ 
(2.7) 

|0). 
(2.8) 



a,p is the creation operator for a boson with on-shell mo- 
mentum p and satisfies the (covariantly normalized) com- 
mutation relations 



= E p S(p-p') 



(2.9) 



For convenience we denote from now on the Lorentz- 
invariant momentum-space integration measure d 3 p/E p 
by dp. 



The assumed factorization in (2.7) of the n-particle 
distribution p(Ci, ■ ■ ■ , Cn) m to a product of single-particle 
phase-space distributions p(Q) for the source current cen- 
ters is often said to represent "independent pion emis- 
sion" In fact, this is only true if the states \n) are 
orthogonal; as pointed out in for the case studied here 
(where J(x) describes a superposition of wave-packets jo 
localized at phase-space points Ci), the emission of a pion 
by one current jo still depends on the positions of the 
other current amplitudes, and if they are closely spaced 
in phase-space, Bose-Einstein symmetrization leads to an 
emission probability which is enhanced by their phase- 
space overlap ("stimulated emission" Q). - The ansatz 
(|2.7| ) also provides a natural starting point for generaliza- 
tions to partially phase-coherent emission, by introducing 
appropriate weight functions into the integrals over the 
phases fa of the source currents pa . 



where Xi is the source center in its rest frame and Aj = 
A(i>i) denotes the Lorentz-boost into the global frame 
with the velocity Vi of the elementary source i. The on- 
shell Fourier transform (p° = E p = \J m? + p 2 ) of J is 
given by 



N 

J(p) = / d 4 xe tp - x J{x) =J2 e 



j (k)= / d^xe^joix) 



<■*■" j^\,-l>>- (2.4) 



As a boosted on-shell momentum the argument of jo in 
( [2.4|) is itself on-shell. We denote the on-shell normaliza- 
tion of j by n : 



d P r i \\ 
~ET \J0(P)\ 



no ■ 



(2.6) 



B. Normalization of the density operator 

The states \n) are not normalized. Their norm 

dp\J{p)\\ (2.10) 



(n \n) = n , n 



depends on the number N of elementary source ampli- 
tudes j and their phase-space positions d and phases 
(2-5) fa (see (2^)). If normalized states were used, this com- 
plicated expression would appear in the denominator of 
(|2.7|) inside the integrals, thereby in general prohibiting 
the analytic evaluation of expectation values 0. The 
only known counterexample which, in spite of construct- 
ing the density operator from normalized states, allows 
for analytical evaluation of the spectra is the coherent 
state formalism ( "covariant current ensemble" , Appendix 
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0) exploited in||o|^|Q. On the other hand, the pre- 
scription (2.7, 2.£) not only permits analytic evaluation 
of the spectra, but also correctly implements stimulated 
emission of bosons, ensuring that for thcrmalized systems 
the Bose-Einstein single-particle distribution is recovered 
in the infinite volume limit 1?|, |25| , ^7|] . 

The normalization of p n therefore mus t be ensured ex- 
plicitly by adding the factor 1/Af(n) in (2.7). The latter 
can be calculated by using 



U J_ J4>i e -i<t>j 



Sij, 



2tt J 2vr 

keeping only the leading terms in N fl20|,[2lf 



(2.11) 



- OO 

■A/(n) ^ 

x f dp 1 ---dp n S„(p 1 ,...,p n ), (2.12) 
S n (Pi,---,P n ) ="o ^PfauPtrJ'-'PiPmPvJ- ( 2 - 13 ) 

er(") 

Here denotes the set of permutations of the indices 
(1, . . . , n), and 

p(p i ,P j ) = — [ d i x 1 fd :i vp(x,v)e~^>-P^ 
n J 

x j * (A{v) ■ Pi ) j (A(v) ■ Pj ) (2.14) 
is the normalized two-particle exchange amplitude: 

dpp(p,p) = l. (2.15) 



The approximation in (2.12) neglects the N n — N(N 



1) ■ • • (N — n + 1) terms in which more than one pion was 
emitted from the same source jo; their relative contribu- 
tion is of order 



N n - N(N -1)---(N -n + l) n(n-l) 



N' 



2N 



(2.16) 



and can be neglected if the mean number of sources is 
much larger than the average number of emitted pairs: 



«JV»»I<n(n-l)) 



(2.17) 



(We denote by (. . .) averages with respect to the pion 
multiplicity distribution p n and by ((...)) averages with 
respect to Pn-) In this limit we thus have 

M{n) = {{N(N - 1) ■ ■ ■ (N - n + 1))) 

X J dp!--- dp n Snip!, . . . ,p n ) . (2.18) 



The limit (2.17) can be enforced by letting in ( |2.q ) 
no -C 1, i.e. each individual source jo emits on aver- 
age much less than 1 pion, and the total emission pro- 
cess is composed of a very large number of elementary 



emission processes. In Sec. IV below we will see that 



the parameter no drops out from the expressions for the 
spectra and correlation functions and thus is physically 
irrelevant. Its sole purpose in the formalism is t o en- 
sure the applicability of the approximations (2.17) and 
( |2.12 ). The corresponding idealization of the emitter as 
a superposition of very many weak source currents should 
be a good approximation for high-energy heavy-ion colli- 
sions, where each event produces anyway a large number 
of hadrons. It may, however, require reconsideration for 
collisions between elementary particles (e.g. e + e~ colli- 
sions at LEP), where only a handful of identical pions is 
created per collision such that the additional smoothcn- 
ing implied by setting n < 1 may be less innocuous. 



C. The ring algebra 

N{n) is related to the "sum of ring integrals" defined 



I 



w(n) = ^| J dp 1 ---dp n S n (p 1 ,...,p n ) , (2.19) 



by 



Af(n) = n\ {(N(N - 1) • • • (iV - n + 1))) w(n) . (2.20) 

The evaluation of uj(n) is done recursively using Pratt's 
"ring algebra" 0. One defines 



Gi(p,q) = n l o J p(P,ki)dkip(ki,k 2 )---dki-ip(ki-i,q) 
G 1 (p,fc)dfc(? i _i(fc, Q ) (2.21) 
with G\ (p, q) = no p{p, q) and normalization 

d = J dpG t {p,p) (2.22) 



(with C\ — no). Identificatio n in (2.21) of p with q 
defines a "ring" ; for this reason (2.22) is sometimes called 
a "ring integral". For the sum of ring integrals u>{n) one 
finds the recursion relation [Hi 



1 - 

;(n) = - Ci uj{n - i) , 

n — ' 



(2.23) 



with uj (I) = no and the definition lu(0) = 1. 

D. The source Wigner density 

The two-particle exchange amplitude p{Pi,Pj) can be 
rewritten in terms of the source Wigner density 17(2;, K): 

p(Pi,Pj) = J d 4 xg(x^( Pl + Pj )) e -<(w-Pi)-* 

d i xg(x,K ij ) e - iq v x . (2.24) 
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Here we introduced the average Kij and the difference 
qij of the two on-shell momenta Pi,Pj. g(x,K) is the 
quantum mechanical analogue of the phase- space distri- 
bution of the emitted pions. According to ( 2.14 ) it is a 
folding of the classical phase-space distribution p of the 
elementarysource centers with their individual Wigncr 
densities |2l{ ]: 

g{x,K)= f d i z 1 i d 3 vp(z,v)g (x-z 1 A(v) ■ K), (2.25) 



go(x,p) = — / d 4 ye-^yf (x + |) j (x - §) . (2.26) 
g and go are normalized to 1: 

J dp d i xg(x 1 p) = J dp d 4 x g (x,p) = 1 . (2.27) 

III. MULTIPARTICLE DISTRIBUTIONS 

The way in which multi-boson symmetrization effects 
manifest themselves in Bose-Einstein interferometry dif- 
fers for ensembles of events with fixed multiplicit y |l| ,|6||8| 
and for multiplicity-averaged ensembles BEl^RjlQ]. In 
this section we provide the relevant formulae for both 
cases. In the multiplicity-averaged case (Sec. pIB| ) we 
allow for arbitrary pion multiplicity distributions; the 
previously studied specific case of Poisson statistics in 
the dilute gas limitEpp,|0,O| is shortly reviewed in 



Sec. IIIC. In Sec. HID we define and differentiate be- 



tween the correlation strength and the incoherence pa- 
rameter and point out that all sources studied here are 
fully chaotic even if the correlation strength differs from 
1. 



It is n orma lized to n(n— 1) ■ • ■ (n—i+l). Note that in the 
limit ( [2.17 ) the dependence on the number of sources, 
N, and its distribution, Pn, drops out completely. This 
differs from the result obtained in ^3 using the so-called 
covariant current ensemble ^0|. The latter is based on 
a density operator which is constructed from a super- 
position of normalized coherent states. In the covari- 
ant current ensemble the one- and two-particle inclusive 
spectra (averaged over the pion multiplicity distribution 
predicted by that model) were found to depend explicitly 
on N and Pjv . This is further discussed in Appendix |A|. 

For the one- and two-particle distributions in an rt-pion 
event Eq. (j3^) reduces to ||| 



_ Jdp2---dp n S n (p 1 ,...,p n ) 
J dpi-~dp n S n (p 1 ,...,p n ) 

^w(n-i) 

..r-\ G »(Pi»Pi)» 



i=l 



u>{n) 



(3.3) 



and 



Nt\ Pl , P2 ) = nin - 1) • Sf---fSJp 1 ,..., Pn ) 

J dpi---dp n S n {p 1 ,...,p n ) 

n , .s i— 1 

%=i v > j=i 

+ <? i (p 1 ,p 2 )G i _ i (p 2 ,p 1 )]. (3.4) 

The 2-particle correlation function from n-pion events 
is given by the ratio 

C<?\ Pl ,p 2 ) = Cf\ q ,K)= ^ )(Pl f£ , (3.5) 

where q = p l — p 2 and K = (p l + p 2 )/2. 



A. Events with fixed pion multiplicity 

We first discuss the case of fixed pion multiplicity n, 
Pm — S mn . With the ingredients from Sec. ||the Lorentz- 
invariant n-pion distribution in an ensemble of events 
containing exactly n pions each can be written as 



Ni n \ Pl ,..., Pn )=E 1 ---E n tx (p n a+ 

n\S n {p 1 ,...,p n ) 



J dpi- - -dp n Snfai, . . . ,p ri 



(3.1) 



It is normalized to n\. The i-particle distri butio n in an 
7i-pion state is (for i < n and in the limit (2.17)) given 
by 



N> n (Pl, ■ ■ ■ ,Pt) = El ■ ■ -EitX (p n o+ • • • 0+0p 4 • 
dp i+1 ---dp n N n n) (p 1 ,...,p n ) 



(n-i)! 

nl Jdpi +1 ---dp n S n ( Pl ,...,p n ) 
(n-i)! Jdpi---dp n S n (p 1 ,...,p n ) 



(3.2) 



B. Multiplicity-averaged ra-particle spectra 

Multiplicity-averaged expressions are now easily de- 
rived by averaging with the pion multiplicity distribu- 
tion p n according to Eq. (2.1). For the n-pion inclusive 
spectrum one finds 



N i (p 1 ,...,p i )=E 1 ---E 



dN 



d 3 pi ■ ■ ■ d 3 p, 
= E 1 ---E i tv (pa+ 



=x;p n jv i (n) ( Pl ,... )Pi ). 



(3.6) 



It is normalized to (n(n— 1) • • ■ (n — i + l)) ( i.e. the i th fac- 
torial moment). Inserting the expressions (|3.3|) and (fTj), 
interchanging the summations over n and i, defining 



Lu(n — i) 
u)(n) 



(3.7) 
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and using the methods of Refs. J^] one finds the follow- 
ing simple expressions: 



dN 



(3.8) 



E\E 2 



dN 



d 3 pi d 3 p 2 



N 2 {Pi,P 2 ) 



= J2 hi +3 [ G i(Pl,Pl)G 3 (p 2 ,P 2 ) 



(3.9) 



+G i (p 1 ,p 2 )G j (p 2 ,p 1 ) . 
For the two-particle correlation function one obtains 



C2(Pl,P 2 ) 



(3.10) 



Si, j h *+j [G l {Pi,Pi)G ] {p 2 ,p 2 )+G l (p 1 ,P 2 ) Gj (p 2 ,Pi)] 
Y^ijhihjGiip^p^Gjip^Pz) 

The corresponding expressions in the previous subsection 
are recovered by setting p m — S mn (i.e. by substituting 



hi = for i > n and hi 



j(n—i) 



for i < n) 



C. "Poisson limit" 

These expressions simplify considerably J3|,|| if the fol- 
lowing special multiplicity distribution is assumed: 



cu(n) N n 



(3.11) 



Here TV is an arbitrary scaling parameter which leaves 
the structure of the theory unchanged. The distribution 
p n depends on v since the sum of ring integrals oj{n) 
depends on v. In the limit of a large phase-space volume 
v (i.e. for a very dilute and sufficiently hot pion gas) the 
higher order ring integrals C(i), i > 1, can be neglected, 
and ui(n) approaches ||||D 



lim uj(n) 



(3.12) 



In this limit p n becomes a Poisson distribution with mean 
multiplicity (n) = uqN: 



w(n) N" 



lim p n = lim - 



-<n> 



(3.13) 



This is why we call this special case the "Poisson limit" , 
with quotation marks to indicate that the multiplicity 
distribution is Poissonian only in the limit v — > 00. In the 
opposite limit p n becomes a Bose-Einstein distribution 
with mean multiplicity (n) = n a N/(l-n N) 0: 



lim p n 

v— *1 



(n) n 



(i + W) 



n+1 



(3.14) 



This can be easily confirmed by using the limiting ex- 
pressions derived in the next section. 

For the choice (3.11) the factors hi in ( p/7) ar e given 
by hi — N 1 and satisfy hi+j — hihj. Eqs. ( 3~8|) ( 3 . 1 C ) 
then reduce to 



N.ip) = H(p,p) , 
N 2 {p l7 p 2 ) = H(p 1 ,p 1 )H(p 2 ,p 2 ) 

+H(p 1 ,p 2 )H(p 2 ,p 1 ) , 



C 2 {p 1 ,P 2 ) = 1 



H(s> x ,Px)H(p 2 ,p 2 ) 
with the auxiliary function 



H(p,q) = J^N i G i (p,q). 



(3.15) 

(3.16) 
(3.17) 



(3.18) 



With the "effective emission function" S(x,K) defined 

by PI 



H iPvP 2 ) 



d 4 xS(x,K)e iq - : 



(3.19) 



the correlation function can then be written in the "stan- 
dard form" [HHUH 



C 2 (q,K) = l 



\Jd 4 xS(x,K) 



fd*xS(x, Pl )fd*yS(x,p 2 ) 



(3.20) 



^i(Pi)iVi(p 2 ) 



Jd 4 xS(x,K) e l "- x 



Jd A xS(x,K) 



Note that this simplification occurs only for this particu- 
lar multiplicity distribution (which in this paper we will 
call "the Poisson limit" for shortness). It allows rela- 
tively easy access to the space-time structure of the ef- 
fective emissi on fu nction, via the last factor in the sec- 
ond term of ( 3.20| ) which can be isolated from the mea- 
sured correlation function by subtracting the 1 and di- 
viding by the (measured) ratio of single particle spectra 
[NiiK)} 2 /[N 1 (p 1 )N 1 (p 2 )}. For other multiplicity distri- 
butions we will see that the space-time interpretation of 
the correlation function is less direct and more compli- 
cated. 

As far as we are aware, except for all previous stud- 
ies of multiparticle symmetrization effects refer either to 
the case of fixed multiplicity or to the above "Poisson 
limit" . 



D. Correlation strength and incoherence parameter 

Since the quantum statistical correlations introduced 
by multiparticle symmetrization effects are expected to 
disappear for large relative momenta, it seems reasonable 
to define the correlation strength X(K) by (!]] 



|<rHoo C 2 (q,K) 



(3.21) 



G 



In the past this correlation strength has often been iden- 
tified with the degree of chaoticity or lack of phase co- 
herence of the particle source: sources which emit par- 
ticles coherently were taken to yield A = while fully 
chaotic sources correspond to A = 1. As we will see be- 
low, this is true only in the "Poisson limit" of the preced- 
ing subsection, but not for general m ultipli city distribu- 
tions where the correlation strength (3.21) can become 
either much larger than 1 or even negative. What re- 
mains true for arbitrary multiplicity distributions is that 
for a coherent source the second (2-particle exchange) 
terms in the square brackets of Eqs. ( |3.4[ ) and van- 
ish, and that for partially coherent sources (like those 
discussed in ]29|-|32]]) they amount to only a fraction of 
the first term even at q — 0. Partially coherent sources 
can thus be characterized by an incoherence parameter 
< Ai nco h(^) < 1 where Ai nco h(-fO is defined by 



Aincoh(-K') 



N° ym (K,K) 



(3.22) 



with A ^ ym denoting the full, symmetrized expressions 
(3.4), Q3.9D, an d N^^™ 1 denoting only the first terms 
inside the square brackets, without the exchange term, 
respectively. 



Since the first and second terms in (3.4) and (3.9) are 
always equal to each other at q = 0, the sources discussed 
in the present paper correspond to 



Aincoh(-f^) = 1 ■ 



(3.23) 



i.e. to fully chaotic sources, for all K and independent of 
the multiplicity distribution. This statement applies also 
to the "Poisson limit" discussed in the preceding sub- 
section, in particular to the limit in which the parame- 
ter no approaches the critical value n c for Bose-Einstein 
condensation as discussed in ||^, [To| , p^ | . Similar to the 
case of an ideal Bose gas in statistical mechanics, Bose- 
Einstein condensation is thus not correlated with the on- 
set of phase coherence in the source. As also pointed out 
m ||], this differs fr om the situation in a laser in which 
the multi-boson state is a single coherent state which ex- 
hibits complete phase coherence, shows no Bose-Einstein 
correlations, and features a Poissonian multiplicity dis- 
tribution. 

Eq. ( [3.2C ) shows t hat i n the "Poisson limit" also the 
correlation strength ( 3.21 ) is X(K) = 1 independent of 
K, and that in this limit it indeed agrees with the inco- 
herence parameter, Ai nco h- 



appeared ||[io|,[l8| which provide a useful basis for com- 
parison. For this reason we adopt the same model here. 

The results of Sees. || and [II show that the multipar- 
ticle spectra are fully determined once the pion multipli- 
city distribution p n has been specified and the functions 
Gi(pi, . . . ,Pj) have been ca lcula ted for all i. From the 
definition of the latter, Eq. (2.21), we see that the Gj are 



completely controlled by the structure of the two-particle 
exchange amplitude p(p, q) or , equ ivalently, of the source 
Wigner densit y g(x , K) (see ( 2.24 )). That the latter can 
(according to (2.25)) be thought of as composed of two 
ingredients, the elementary source currents jo and the 
phase-space distribution p(£) of their centers, does not 
matter in practice, i.e. it has no observable consequences 
Fo r a given source Wigner density g(x, K), equa- 



tion (2.25) in general allows for infinitely many different 
decompositions into a classical phase-space distribution 
p and an elementary source Wigner density go- This be- 
comes manifest in the case of nonrelativistic Gaussians 
(studied e.g. in ]7|,p|,^3| [37| ) where the Gaussian widths 
of g(x, K) in coordinate and momentum space, R and 
A, are related to those characterizing p (Rq and Ao) and 
go (a/y/2 and l/(v / 2c) — in those papers go was con- 
structed from minimum uncertainty wave packets of size 
a) by 



E>2 d2 

It — _Kq 



a 

Y 



A 2 



2a 2 



(4.1) 



Obviously this has, for given R and A, an infinity of 
solutions (i?0j Ao,ct). Therefore, the wave packet width 
a is not directly measurable and can at most be bounded 
by the two measurable quantities R and A. 



A. Static sources 

We will here study the non-relativistic Gaussian source 
Wigner density 



l(x,p) 



E p 6(x°) 
(2ttRAY 



exp 



2A 2 



(4.2) 



which describes instantaneous emission of particles with 
a Gaussian (non-relativistic thermal) momentum distri- 
bution from a sphere with a Gaussian density distribu- 
tion. The momenta do not depend o n po sition, i.e. x and 
p are uncorrelated. For the source ((ij) the two-particle 
exchange amplitude reads 



IV. GAUSSIAN SOURCES 

Due to their analytic tractability Gaussian source mo- 
dels have always been popular in Bose-Einstein interfero- 
metry (for recent reviews see p2|). It was recently dis- 
covered [Q that they also allow for an analytic evaluation 
of multi-particle symmetrization effects. Since then se- 
veral investigations exploiting this analytic solution have 



P(P1,P2> 



V Ei E2 
(27rA 2 )i 



exp 



K 

~2A 2 



R 2 q 2 



(4.3) 



The higher-order exchange amplitudes G n (p 1 ,p 2 ), n>l, 
can be calculated ana lytica lly from recursion relations 
§ generated by Eq. ( ^2l|) : 
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Tl'rt c 



o ^ n 



exp 



Gn{Pl,P 2 ) 

VWE^ ~ (2ttA 2 ) 
Rl = a n R 2 , A 2 n = a n A 2 

__ 1 (v+ l) n + (v- l) n 
" _ v (v + 1)™ -(v- 1)" 

2 2n v 



K 

' 2A 2 



lt„ q 



< 1, 

3/2 



+ 1) 2 » - (« - l) 2 



< 1 . 



where 



2.RA > 1 



(4.4a) 
(4.4b) 
(4.4c) 

(4.4d) 
(4.5) 



is a measure for the total phase-space volume occupied by 
the source Wigner function g(x,K). (The latter would 
be given by v 3 , and v=l corresponds to the minimum 
value allowed by the uncertainty relation.) The ring in- 
tegrals C„ are given by 



C n — C n dn 



(4.6) 



It is easy to generalize expressions (4.2)— (4.5) to the case 
of a Gaussian source with different widths Ri,R 2 , R3 and 
A 1 ,A 2 ,A 3 in the three spatial directions, with corre- 
sponding phace-space volume factors Ui, i>2, i>3. 

From the above results one sees that the parameter 
no (the normalization (|]^) of the elementary sources) 
drop s out from all observables: Eqs. (4.6), ( p. 23 ) and 



G 



imply that u>{n) 



and hi 



Since also 

1, al l factors of no cancel in Eqs. (3.3), (3.4), and 



( 3.8)— (3. 1C ) . We will therefore, without loss of generality, 
set n = 1 from now on J38j. 

In Eqs. (fhj) we have rewritten the results of Zimanyi 
and Csorgo |7j in such a way that the essential structure 
of the multipartic le sym metrization effects becomes ma- 
nifest: first, Eq. (4.4a) shows that n-body symmetriza- 
tion effects contribute to the multiparticle spectra with 
relative weight c„<l, and these contributions correspond 
to effective source functions with modified widths R n <R 
and A„<A. This reflects the fact that bosons like to 
cluster in phase-space. Second, the weights c n and the 
strength of the clustering, given by controlled by 

a single parameter, the available phase-space volume v. 

For large v multiparticle effects are weak; one finds for 
n> 1 



lim a n 

v — >oo 

lim c n 

v — >oo 



1 

3 

71 

712 \V 



lim C n = — o ( - 



3(n-l) 



3(n-l) 



lim coin) 



1 



(4.7a) 
(4.7b) 

(4.7c) 
(4.7d) 



In this limit the weight of n th -order symmetrization ef- 
fects is suppressed by the n— 1 st power of the input phase- 
space density do = l/v 3 (i.e. the phase-space density in 
the absence of multi-boson effects). 



As v becomes small and approaches the uncertainty 
limit, v—>-l+e, one finds for n > 2 



lim a n = 1 — e, lim " +1 = 1 — 2 ( — 
lim c„ = 1 — |(n— l)e , 

u->l+e 2 

lim C n = 1 — Ine , 



l'-»l + e 

lim 



lim u)(n) = lim c n — 1 ~ |(n— l)e. 



v— >l+e 



(4.8a) 
(4.8b) 
(4.8c) 
(4.8d) 



(Of course, ai=Ci=C%=l independent of v.) Now all n- 
boson terms contribute with equal weight. However, the 
bosonic clustering effects become weak, a n — >1, since the 
source is already as small as allowed by the uncertainty 
principle; all n-boson terms thus reflect the same effective 
source parameters. In the following section these features 
and their implications for Bose-Einstein interferometry 
will be studied quantitatively. 



B. Gaussian sources with a>p-correlations 

In this subsection we shortly describe a Gaussian 
model, first suggested by Zajc in which the momenta 
of the emitted particles are correlated with the emission 
points, as is the case e.g. in hydrodynamically expand- 
ing sources. The model allows to control the strength 
of the x-p-correlations by tuning a parameter s and thus 
presents a valuable playground to investigate the quali- 
tative effects of £-p-correlations on HBT interferometry. 
We will show that, for the purposes of studying multi- 
boson effects, it can be completely mapped onto the static 
Gaussian source presented in the previous subsection: af- 
ter a suitable redefinition of the physical parameters, the 
algebra for calculating multi-boson effects becomes iden- 
tical to the already solved static problem. 

Zajc's model starts from the non-relativistic classical 
phase-space distribution 



p(x,p) 



where 



R s 



S(x°) 



(27Ti? s A ) 3 



x exp 



(4.9) 



2(l-s 2 ) 



2s 



x ■ p 

Ro&o 



El 

A 2 



Rq\/1-s 2 , < s < 1 and p = 



mv. 



(4.10) 



The term proportional to s in the exponent couples mo- 
mentum to position; in the limit s — * 1 one obtains per- 
fect x-p-correlations: 



lim p(x,p) ~ S(x ) 5 

s— *1 



Ro 



P 

A 



(4.11) 



This violates, of course, the uncertainty relation; the 
threshold at which the effective phase-space volume fac- 
tor v s = 2i? s A becomes smaller than 1 is given by [p7| 



8 



'crit 



1 - (2i? A )" 



(4.12) 



To obtain a valid source Wigner density which does 
not violate the uncertainty relation, we now fold |21j] p 
according to ( 2.25 ) (leaving out the factor 7 4 in the inte- 
gration measure since the source is non-relativistic) with 
the Wigner density of the elementary source currents jo, 
which are taken as non-relativistic minimum uncertainty 
Gaussian wave-packets of width a: 



9o(x,p) 



E p S(x°) 



exp \_-x 2 jo 2 - a 2 p 2 ] 



(4.13) 



After lengthy, but straightforward algebra we thus obtain 
the source Wigner density 



i \ E p 
9{x,p) = 



d(x°) 



(27Ti?A) 3 



exp 



(x — sL 2 p) 2 p 2 



2B 2 



2A 2 



(4.14) 



where [B9fl4(J 



L 2 = 



RqAq 
A 2 



R 2 A 2 = R 2 A 2 - s 2 R 2 A 2 



1 

> - . 
~ 4 



0~ 

Y 



A 



1 

2^2 



(4.15a) 
(4.15b) 
(4.15c) 



Eq. ( 4.14 ) differs from ( |4.2| ) only by a p-dependent shift 
of the spatial coordinate x. The new length scale L char- 
acterizes this x-p-correlation. Otherwise the two expres- 
sions are completely identical, aft er re defining the physi- 
cal parameters R, A according to ( 4.15 ). The correspond- 
ing two-particle exchange amplitude is given by 

'2 d2„2 



(27rA 2 )i 



exp 



K 

2 A 2 



- isL 2 K-q 



(4.16) 



This differs from (4.3) only by the phase factor inv olving 
L 2 . Since K q = {p\ —p 2 )/2, in the ring integrals (2.21) 
the phase factors involving the integration momenta can- 
cel, leaving only those corresponding to the external mo- 
menta. Since the particle spectra are real and always in- 
volve products of such pha se fact ors with their complex 
conjugates (c/. Eqs. ( |Q| ) and (|3.9|)), they completely 
drop out from the calculation of higher-order multi-boson 
effects. The algebra of the previous subsection thus goe s 
through without modification; the s-dependence of (4.S) 
which describes the strength of x-p-correlations in the 
source has been completely absorbed into the redefined 
source parameters R, A. 



V. RESULTS 

In this section we give a comprehensive and quantita- 
tive discussion of the behaviour of the correlation func- 
tions, the correlation strength, and the HBT radius pa- 
rameters, both for events with fixed multiplicity and for 



multiplicity-averaged event ensembles. The results for 
fixed multiplicity can be obtained in two different ways: 
either by using Eqs. (^)-(3.5), or from the expression 
(|3.10D for multiplicity -averaged events with the special 
multiplicity distribution p m =S rnn . We have checked that 
both paths lead to identical results, but the second path 
allows for an easier presentation. 



A. Correlation functions: Bose-Einstein and residual 
correlations 



We begin by rewriting ( 3.1 0| ) as follows: 



C 2 (Pi,P 2 ) 



Ei,,=i h i h j GiiPnPi) Gj(p 2 ,P 2 ) 



' x Sjj=i h i+j G i (p 1 ,p 2 ) G j (p 2 ,p 1 ) 



ct S \q,K)[l + R 2 (q,K) 



(5.1) 



The second term in brackets, R 2 , arises from Bose-Ein- 
stein symmetrization of the two measured momenta, p 1 
and p 2 . As can be seen from the discussion in Sec. Ill C| , 
in particular Eq. ( |3.19 ), this is the term which contains, 
via a Fourier transform, information about the effective 
source size; to extract this information is the goal of HBT 
interferometry. The expression inside the bracket shows 
manifestly that th e sour ce is fully chaotic, as discussed 
at the end of Sec. HID : at q=0 we have R 2 = 1, and 
both terms inside the bracket beco me equal. 

The interesting feature of Eq. (5.1) is the prefactor 



C 2 es \q, K) which, following Zajc fllLM who first noticed 



this effect, we call "residual correlation". Equation (5.1) 
clearly exhibits the physical origin of the residual corre- 
lations: they describe the result of the different combina- 
torics involved in creating background pairs from mixed 
events (the denominator of C^ 10 ^) and "sibling pairs" 

from the same event (the numerator of C^ 6 ^), without 
symmetrizing with respect to the two measured momenta 
p 1 and p 2 . [Note that the measured momenta p Y and p 2 
are still symmetrized with respect to all other momenta 
of particles within the respective events; for this reason 
C^ 08 ^ can, unfortunately, not simply be interpreted (and 
measured separately) as the ratio of unlike-sign pairs 
from the same event divided by the like- or unlike-sign 
background pairs from mixed events: in that ratio the 
numerator involves symmetrization only with respect to 
the momenta of those particles having the same charge 
as the measured ones, leading to different combinatorics.] 
For sources of the type described in Sec. [n| the residual 
correlations are absent if and only if the multiplicity dis- 
tribution corresponds to the "Poisson limit" of Sec. [Ill C 
(see Eq. ( 3.20| )) where hi+j=hihj. And only in this limit 
it is possible to isolate the information on the source ge- 
ometry by extracting from the Bose-Einstein correlator 







i?2 a simp le r atio of particle spectra (see again (3.20)). 
Equation (5.1) thus shows that the usual formalism for 
HBT interferometry must be significantly modified if one 
allows for arbitrary multiplicity distributions in a situa- 
tion where multi-boson symmetrization effects must be 
taken into account. This is true in particular for sources 
with a fixed pion multiplicity. While this need for a mod- 
ified extraction procedure for the HBT radii was noticed 
and a specific suggestion for dilute sources was made in 
Jl3|, we present here a completely general approach to 
this problem. 

The goal of the remainder of this section is to give 
a quantitative assessment of the importance of residual 
correlations and multi-boson effects on HBT radii, by ex- 
ploring the full physical parameter range for the Gaussian 
sources discussed in Sec. ^ and for a variety of standard 
multiplicity distributions. This study reveals a number 
of unexpected new features of two-particle interferometry 
which have not been previously investigated. As we will 
show, they can be used in new ways to obtain information 
on the phase-space density of the source at freeze-out. 
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FIG. 1 
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Correlation function Ci (solid), residual cor- 
(dotted), and Bose-Einstein correlation R2 
(dashed) for events of fixed multiplicity n — 10 (top) and 
n = 50 (bottom), for vanishing pair momentum k. (Note 
that for k — the correlation function is isotropic in q.) In 
both plots the phase-space volume factor is set to v = 2.56. 

It is convenient to introduce scaled momenta, 

k = K / A , q = Rq, (5.2) 



and write k = \k\, q — \q\. Pair momenta are thus mea- 
sured in units of the width A of the source in momentum 
space, while relative momenta are scaled by its width R 
in coordinate space. 

Explicit expressions for the correlation function (pd]) 
for a Gaussian source are given in Appendix |b| Fig. |l| 
shows two sample correlation functions for events with 
fixed multiplicities n = 10 and n = 50, respectively. Ob- 
viously, in both cases the correlation strength A (defined 
as the ratio of the correlation function at q = and 
q = 00 minus 1) is less than one. For the selected val- 
ues of the pair momentum k and the phase-space volume 
1?, the correlation strength actually vanishes at n = 50. 
Fig. H shows that this is caused by the residual correla- 
tions which become stronger with increasing phase-space 
density n/v 3 , pulling the intercept at q = down. This 
effect is studied in more detail in Sec. |V Q 

For the case of fixed multiplicity studied in Figures |l| 
the residual correlations are negative, i.e. the correlation 
function approaches its asymptotic value at q — ► 00 from 
below. This was found before in Fig. 6 of the last paper of 
Refs. H (although not explicitly noted in the text), and 
the origin of the effect was correctly identified in |L3| and 
independently in p5fl . We want to emphasize, however, 
that the negative sign of the res idua l correlations is not 
generic: as will be seen in Sec. V C , it depends on the 
multiplicity distribution of the emitted particles, and a 
Bose-Einstein distribution, for example, yields positive 
residual correlations. 



B. Multiplicity distributions 

Before proceeding further, we specify in this subsection 
a few particular forms for the measured pion multiplicity 
distribution. They will be used below to illustrate our 
general results with specific examples. In particular, we 
will show results for: 



Pn 
Pn 



(n) n 



(fixed multiplicity m) 



-<n) 



p(n,k) 



(Poisson distribution) 
(n+fc-l)! 



n\ (fc-1)! 



(negative binomial distribution) 



pr{n,k) = 



(l+(n)T +1 ' 

(Bose-Einstein distribution) 
1 f kn ' ' 



,-kn/(n) 



nT(k) 
(Gamma distribution). 



(5.3a) 
(5.3b) 

(5.3c) 
(5.3d) 
(5.3e) 
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The Poisson distribution (5.3b) is ob taine d from the 
negative binomial (NB) distribution ( 5.3c ) by setting 
k — > oo; the Bose-Einstein distribution (5.3c) is obtained 
for k—1. For the Poisson distribution one has 



(n(n-l)) = (nf 



(Poisson) 



(5.4) 



while the Bose-Einstein distribution gives a larger vari- 
ance: 

(n(n-l)) = 2 (n) 2 <*=^ (n 2 ) - (n) 2 = (n) (l + (n)) 
(Bose-Einstein) . (5-5) 

For the case ( |5.3a ) of fixed multiplicity the variance van- 
ishes, of course. 



C. Normalization and correlation strength 

The Bose-Einstein correlation R 2 equals 1 at q — and 
vanishes at q — > 00. The asymptotic normalization of the 
correlation function C2 and the value of the correlation 
strength are thus controlled by the residual correlations 

£f(res) 

The normalization Af does not depend on the pair mo- 
mentum K] the results of Appendix [b] give 

h 2 

. MC/.K) = — 

q — >oo 



Af = lim C 2 (q,K) - 

Oi) 



(5.6a) 



lim N 

v — >oo 



(n(n— 1)) 



2 Bose-Einstein . 
1 Poisson , 



(5.6b) 



lim AT = 1 P ° ^ 



1 Bose-Einstein , 

l-«n) + l)e-<"> D „._ (5.6c) 
(l-e-<">) 2 



Poisson 

For fixed multiplicity n this reduces to 

^•(n) = u{n)u>(n-2) 
u> 2 (n — 1) 



lim A/" ( 



lim AA (n) = 1 



n) _ n{n- 1) 



(5.7a) 

(5.7b) 
(5.7c) 



The generic difference of Af from unity can be clearly 
seen in the upper panel of Figure [l]. 

For the correlation strength one finds, using results 
from Appendix [b|, 

(hi) 2 A(k) 1 



A(k) = 2 ■ 
lim A(k) = 1 



ha B(k) 



limA(K)= (1 ' P ° )2 ^"r^-l 

f 1 



,-{w) (n)-l+e-<"> 
l-«n) + l)e-<"> 



(5.8a) 
(5.8b) 

(5.8c) 

Bose-Einstein , 
Poisson . 



For fixed multiplicity n this reduces to 

lim AW((t) = 1 > lim A (n) («) = --■ (5.9) 

v— >oo 11— >1 n 

One sees that for fixed n in the limit v — > 1 (i.e. for 
small sources which saturate the uncertainty relation) the 
residual correlations become so strong that the effective 
correlation strength A becomes negative. 




v=2RA 

FIG. 2. Correlation strength X^(k) as a function of the 
dimensionless phase-space volume v — 2RA/H. Results are 
shown for two fixed multiplicity values (n = 50 (solid) and 100 
(dashed)) and for three values of the reduced pair momentum 
K = K/A = 0, 1, and 5 (from right to left). 



Closer inspection of Eq. ( |5.8c ) for the Poisson multi- 
plicity distribution shows that it is not really the phase- 
space volume factor v, but the average phase-space den- 
sity 



(n) 

1 1 3 



(5.10) 



which controls the strength of multi-boson effects. For a 
source of minimal phase-space volume v=l, the correla- 
tion strength A approaches zero as (n)—>oo (large phase- 
space density d — > 00); this agrees with the expectation 
of strong multi-boson effects in this limit. On the other 
hand, in the opposite limit (n) — > (small phase-space 
density d — * 0) A approaches unity, the same value as 
for sources with large phase-space volumes v — > 00. For 
strong multi-boson effects to arise, it is thus not suffi- 
cient that the phase-space volume factor v is small, but 
the phase-space density d = (n) / v 3 must be large. 

Let us now discuss some numerical results. In Figures || 
and H we show the correlation strength A as a function of 
the phase-space volu me f acto r v = 2RA/H. Whereas, ac- 
cording to Eqs. (5.8) and J5.9| ), the correlation strength is 
independent of the scaled pair momentum k in the limits 
v— >oc and v— >1, Figure || shows that in the intermediate 
range 1 < v < 00 the correlation strength A'™) develops 
a K-dependence. At fixed phase-space volume v, the ef- 
fective strength of the correlations diminishes for smaller 
k and/or larger multiplicity n. 

The dependence of the correlation strength on the mul- 
tiplicity distribution is shown in Figure ||. For fixed 
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k = 0.3 and constant mean multiplicity (n) — 18, Fig. || 
gives A as a function of v for four different multiplicity 
distributions. One sees that for different multiplicity dis- 
tributions the residual correlations manifest themselves 
quite differently. There is no universal scaling with the 
average phase-space density d defined above; however, 
for typical momenta k ~ 1, significant residual corre- 
lation effects (i.e. noticeable deviations of A from its 
asymptotic value lim^oo A=l) are seen for phase-space 
densities above 0.2-0.3 particles per unit phase-space vol- 
ume, d > 0.2 — 0.3. At lower momenta k < 1 similarly 
strong effects arise from even smaller values of d, whereas 
only very high phase-space densities lead to strong multi- 
boson effects at large pair momenta n ^> 1. These con- 
clusions agree with those of Lednicky et al. Il3| . 
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FIG. 3. Correlation strength A(k) as a function of the di- 
mensionless phase-space volume v = 2RA/H, for fixed pair 
momentum k=0.3. Results are shown for four different pion 
multiplicity distributions (Poisson, Bose-Einstein, Gamma 
distribution with fc=8, and negative binomial (NB) distribu- 
tion with fc=6). In all cases we fixed (n)=18. 

Note that the multiplicity distributions studied in 
Fig. U refer to the observed pion multiplicity, after full in- 
clusion of multi-boson effects; the measured multiplicity 
distribution is kept fixed while varying the phase-space 
density. The presentation in Fig. |3| differs from the dis- 
cussion of the "Poisson limit" in Sec. [Ill C where the in- 
put distribution was fixed to be Poissonian, an d the mea- 
sured multiplicity distribution p n {v) of Eq. (3.11) then 
depended on v. The distribution p n {v) corresponds to 
a constant correlation strength A=l and would in Fig. || 
interpolate as a constant horizontal line between the Pois- 
son distribution at v— >oo and the Bose-Einstein distribu- 
tion at v— >1. This must be taken into account when com- 
paring with results from the "Poisson limit" like those 
discussed in ff^|lO|JiT 



D. HBT radii 

The information on the space-time structure of the 
source is contained in the Bose-Einstein correlation func- 
tion R 2 (q,K) of (5.1). The discussion of the "Poisson 



limit" in Sec. [HIC| (cf. Eq. ( |3.20| )) showed that this in- 
formation is "contaminated" : the two-particle exchange 
term ^ hi + j Gi(p 1 ,p 2 ) Gj(p 2 ,p 1 ) is normalized by the 
corresponding "direct term" evaluated at the observed 
momenta p 1: p 2 (i.e. by £ h i+j G l (p 1 ,p 1 ) G j (p 2 ,p 2 )) 
instead of the same expression evaluated at the average 
momentum If (i.e. J2 hi+j Gi(K, K) Gj(K, K)). Only 
in the latter case would R 2 {q 1 K) allow for a "clean" ex- 
traction of the space-time structure of the source. For 
heavy-ion collisions where the sources are large it is 
known jl2|,[l3| that (at least in the "Poisson limit" ) these 
contamination effects from inadequate normalization are 
negligible. For small sources, like those created in e + e~ 
collisions, they may, however, seriously affect the extrac- 
tion of the source radius. 



In the "Poisson limit" Eq. (3.20) shows that this type 
of contamination can be avoided by first extracting from 
R 2 a (directly measurable) ratio of single part icle spectra. 
For a dilute Gaussian source of the type (fh|) (i.e. in the 
limit v — > oo) this ratio is given by 



[Ni(K)} 



Ni(Pi)Ni(p 2 ) 



exp 



^4A 2 



= '" MM ^ 



(5.11) 



For arbitrary multiplicity distributions the residual corre- 
lations prohibit such a simple procedure (the denomina- 
tor of i?2 cannot in general be written as a product of sin- 
gle particle spectra), and the contamination of the space- 
time information about the source cannot be avoided. 
The results presented in this subsection allow, however, 
to estimate the relative magnitude of this contamination 
effect. 

As usual, we parametrize the correlation function 
R 2 (q,K) by a Gaussian in q with JC-dependent width 
parameters (HBT radii). The spherical symmetry of the 
Gaussian source ( |4.2| ) ensures that the HBT radii de- 
pend only on the modulus n of the scaled pair momen- 
tum (see Appendix ^). However, since the denomina- 
tor o f R 2 invol ves the two measured momenta p 1 and p 2 
(see (5.1) and (|B4|)), the correlation function R 2 depends 
separately on the components q\\ and q± of q which are 
parallel and orthogonal to K . This is an artifact of the 
above-mentioned "contamination effect" ; the difference 
of the width parameters in the q\\ and q± directions thus 
is a quantitative measure for this undesired contamina- 
tion. 



We thus write 

R 2 (q,K) =cxp [-q\ R\{K)-q\ R\(K) 
q±=qsm6, q\\ = q cos6 , 9 = L{q,K) 



(5.12) 
(5.13) 



The HBT radii are then obtained by separating the con- 
stant and ^-dependent terms in 



dq 



1^2 Mq,k) 



= Rl(K) 



(5.14) 



9=0 



>s 2 6 (r\{K) - R 2 ± (K) S 
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The correspond ing explic it e xpressions are given in Ap- 
pendix^, Eqs. (CI) and (p2|). Before proceeding further 
it is useful to establish a baseline in the absence of multi- 



boson effects, by keeping in (CI) and (C2) only the lowest 
order terms i=j=l. Denoting these reference values by 
a tilde we find 



Rj_ = ij| = B* ( l 



1 



(5.15) 



independent of the pair momentum k. The HBT radii are 
seen to be smaller than the Gaussian input radius R by a 
correction ~ 1/w 2 ; this is the "contamination" from the 



ratio (5-11) of single particle spectra mentioned above. 



In the following we will discuss multi-boson effects on 
the HBT radii (including additional multi-boson-induced 
"cont amination" effects) as deviations from the baseline 
(PHI) . To this end we define the dimcnsionless ratios 



n 



Rl 



R 2 1 



ti n Xtj^ 

R 2 



(5.16) 



Like the HBT radii themselves, these ratios depend on 
v, k and (n) via multi-boson symmetrization effects; in 
the absence of such effects (in particular in the limit v — > 
oo) they approach r\ = 1 and r 2 = 0. In particular, 
r 2 provides a relative measure for multi-boson-induced 
spectral contamination effects. 




FIG. 4. The ratios n (thick lines) and T2 (thin lines) as 
functions of v = 2RA/h, for fixed multiplicity n = 20 and 
for three values of the scaled pair momentum k = 0, 1, 5. ri 
vanishes identically for n = (see Eq. (jC2|) ) . 

Analytic expressions for r% and r 2 in the quantum sat- 
uration limit v — ► 1 (in which the pions in the source all 
occupy the same quant um state and form a Bose conden- 
sate) are given in Eqs. (|C3|) . Note that in this limit both 
Rj_ and Ru vanish even in the absence of multi-boson ef- 
fects, due to the factor (1 — \) in ( 5.15 ). This leads to 
a flat Bose-Einstein correlation function 1 + R 2 . In 
this flatness of the correlation function was incorrectly 
attributed to the onset of phase coherence, resulting in 
an absence of Bose-Einstein correlations at low values of 
q. Lednicky et al. p3[ clarified this error by pointing 
out that even in this limit the Bosc-Einstcin correlations 



are still present (the exchange term R^ = 1 is as large 
as the direct term), and that the misunderstanding in 
resulted from an incomplete analysis of the normal- 
ization of the correlation function in the Bose-Einstein 
condensation limit. The latter was also discussed in [n0|. 
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FIG. 5. Top: The ratio ri as a function of the scaled pair 
momentum k, for two fixed multiplicities (n = 10 (thick lines) 
and n = 20 (thin lines)) and three values for the phase-space 
volume factor v = 1,2,5. For v = 1, r\ is independent of k 
and less than 1; but for v = 1 + e, n will eventually again 
approach 1 at k — » oo. Bottom: The same for the ratio T2\ it 
vanishes for v = 1. 



Eqs. (C3) show that for all considered multiplicity dis- 
tributions lim.y_>i ri < 1. The numerical results below 
show that the same holds true for all values of v. In other 
words, multi-boson effects always reduce the HBT radii: 
R\ < R 2 (l - 4j). For v > 1 this is largely the conse- 
que nce o f Bosc-clustering in the source — QnR^ — 
see (4.4)) and confirms similar findings in many earlier 
papers. However, for v = 1 the source is already as small 
as allowed by the uncertainty relation and thus cannot be 
compressed any further by Bose-clustering. In this case 
a reduced value n < 1 must be blamed entirely on multi- 
boson contributions to the spectral contamination effect 
from the product of single particle spectra in the denomi- 
nator of the Bose-Einstein correlator R 2 . Eqs. (C3) show 
that these contributions disappear for vanishing phase- 
space density d = (n}/v 3 — ► 0, even in the limit v — ► 1: 



lim lim r± = 1 



(5.17) 
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(For fixed multiplicity the lowest allowed value is n — 2 
for which ri satisfies the same rule.) However, for large 
phase-space densities d — > oo they are large and the ratio 
T\ vanishes: 



lim lim T\ 

(n>— >oc v— >i 



lim lim 



(n)-oo«-*l i? 2 (1 - 4y) 



0. 



(5.18) 
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FIG. 6. Top: The ratio ri as a function of v, for vanishing 
pair momentum and fixed average multiplicity (n) — 20, for 
four different multiplicity distributions. Bottom: The same 
for the ratio r%, but for k=1 (since r?=0 for k=0). 

These last results can be generalized to arbitrary values 
of v by rephrasing them in terms of d: 



lim n = 1 , 
lim T2 = , 
lim ri = lim T2 = . 

d — >oc d — >oo 



(5.19a) 
(5.19b) 
(5.19c) 



For the limit d — > oo the proof relies on the fact that, 
at fixed i;, this requires (n) — > oo; therefore only the 
asymptotic behavior of (4.4) for large n is needed, in 
particular 



lim a n 



(5.20) 



(Incidentally, this implies linin^oo R n A n = ^j- = ^: 
in the limit n — oo the effective source corresponding 
to GniVxzJh) saturates the uncertainty relation.) Equa- 
tions ( 5.1 9c[ ) then follow easily from ( |Cl| ) and (|C2|). 



Figures ^-|^ show numerical results for r\ and T2 for 
a variety of multiplicity distributions as functions of the 
phase-space volume factor v and the scaled pair momen- 
tum k. In Figures ^ and ^| we study events with fixed 
multiplicity (n = 10 and n = 20). Multi-boson sym- 
metrization effects on the HBT radius R± and their con- 
tributions to the difference — R\ (which is entirely 
due to the above-mentioned "spectral contamination ef- 
fect") are seen to increase as the phase-space density d 
increases and the scaled pair momentum k decreases. On 
the other hand, Fig. |^ shows that for large pair momenta 
k the multi-boson effects vanish, irrespective of the value 
of v. As noted in |ll| this reflects the exponential de- 
crease of the momentum-space density ~ exp(— k 2 /2) at 
large particle momenta. This also explains why increas- 
ing k compresses the curves in Fig. ^ to the left. 
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FIG. 7. Top: The ratio n as a function of k, for a fixed 
phase-space volume factor v = 5 and fixed average multi- 
plicity (n) = 20, for four different multiplicity distributions. 
Bottom: The same for the ratio ri. 

Note that the multi-boson effects set in rather sud- 
denly. (As discussed in Appendix the vanishing of T2 
in the limits v — ► 1 and/or k — > is due to its definition 
and does not imply the disappearance of multi-boson ef- 
fects in these limits.) Unfortunately, we have not been 
able to find a simple scaling law for the threshold at which 
multi-boson effects on the HBT radii set in and which 
properly describes the observed qualitative trends. 

In Figures || and we study the dependence of 
multi-boson symmetrization effects on the measured pion 
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multiplicity distribution. While the dependences on 
v (Fig. ^) and on k (Fig. are qualitatively simi- 
lar to the case of fixed multiplicity on a quantitative 
level one observes quite significant variations for differ- 
ent multiplicity distributions at fixed mean multiplicity 
(n). This demonstrates that there is no simple scaling 
of multi-boson effects with the average phase-space den- 
sity d = (n) /v 3 . Of the four distributions considered 
here, the Bose-Einstein multiplicity distribution gives by 
far the strongest multi-particle symmetrization effects on 
the HBT radii: it is less effective in suppressing high- 
multiplicity contributions than the other considered dis- 
tributions. 



E. The range of the residual correlations 

Looking at Figure [l] and Appendix [b| one sees that 
the residual correlation C^ ^ (q, K) can also be para- 



metrized by a Gaussian: 



A(k)-1 p-«l r l( K )-Q?i r n(«) 
2 ° 



(5.21) 



The two newly introduced radius parameters t±(k), ru (n) 
can be interpreted as the "range" of the residual correla- 
tions in coordinate space. A Gaussian parametrization of 



the full correla tion f unctio n C^j q, K) in (5.1) is obtained 
by combining ( 5.21 ) with ( 5.12 ) . 

From the ansatz (5.21) one obtains the "residual corre- 
lation radii" by separating the constant and ^-dependent 
terms in 



R 2 



R 2 



B(k) 



AfB{ K )-A( K ) dq 



r,2 



(5.22) 



9=0 



where we used ( |5.6a| ) and ( 5.8a ) to rewrite the coefficient 
in front of the g 2 -derivative. T he e x plici t expressions 
are given in Appendix |c|, Eqs. (|C^)-( |C9|) . Using the 

and 



asymptotic expressions (p9h we find from 




lim ri(/t) = lim rjf(«) = 0; (5.23) 

V — >00 V — >OG 

for dilute sources the residual correlations disappear, as 
expected. In the opposite limit v — > 1 the two "residual 
correlation radii" are nonzero, but again equal to each 
other: 



lim(ri(«)-rf(K)) =0. 



(5.24) 



What matters in practice is the ratio of the "residual 
correlation radii" to the HBT radii, r±/R±. For fixed 
multiplicity n, this ratio can be calculated analytically 
in the limit v — > 1: 



lim 



ri(«) 



1 n- 1 



R]_(k) 4n+l 



(fixed n). 



(5.25) 



For n — > oo this approaches the value \ from below; for 
intermediate values 1 < v < oo the ratio is between 
and j. The same is true for the Poisson multiplicity 
distribution for which (C9) gives in the limit of large (n) 



lim 



.(«) 



1 

4 



-2 



\n) 



1 



+ 0(e- {n) ) (Poisson). (5.26) 



From these analytical results it seems that the "range" 
of the residual correlations in coordinate space is always 
less than half the effective source size as measured by the 
HBT radii. Below we confirm this numerically. Thus, 
with two well-separated length scales involved, strong 
multi-boson effects should be clearly recognizable in the 
2-particle correlation function, as exemplified in Fig. |l| 
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FIG. 8. Top: w-dependence of the ratio rj_/R± at vanish- 
ing pair momentum k = 0, for a variety of multiplicity dis- 
tributions, keeping the average multiplicity fixed at (n) = 20. 
Bottom: The same for {r^ — r 2 ± )/R? at k= 1. The grey dashed 
lines indicate an interpolation as described in the text. 



In Fig. U we show the ratios r\jR\ and (r 2 — r\)lR 2 
as functions of v, for the same values of k and (n) as 
m Fig. |. One sees that indeed always r± < and 
that the ratio of these two parameters is largest at small 
values of v. The ratio (r 2 — r'jJ/R 2 alternates sign, but 
its absolute value alw ays st ays below about 5% and also 
below the ratio ri of ( 5.16 ) (see Fig. ||). 
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We note that the range parameters r\{n) and rjj(«) 
become ill-defined when the c orrelation strength A(k) 
passes through 1 (see Eq. (5.21) and Fig. ||). We checked 
that at this point the function C™ s (q) is approximately, 
but not quite exactly constant. The remaining, albeit 
very weak g- depe ndence cannot be properly absorbed by 
the Ansatz (5.21), resulting in a divergence of the ex- 
pression (5.22). Since this is an artefact of the Ansatz 
(5.21) and does not reflect any really singular behaviour 
of the residual correlations, we smoothed this singularity 
in Fig. |^ by interpolating by hand (grey dashed lines) 
across the point X(k) = 1. 



VI. EXPERIMENTAL CONSEQUENCES 

Figure [l] shows that, for fixed n, the residual corre- 
lations reduce the correlation function while the Bose- 
Einstein correlations increase it. The interplay of these 
effects leads to a minimum of the correlation function at 
intermediate values of q. This is different from Coulomb 
or strong final state interactions (FSI) which are concen- 
trated much more strongly near q = (and which in prac- 
tice one can correct for p2[). A minimum at nonzero q, 
after FSI correction, thus points to the presence of resi- 
dual correlations arising from multi-boson symmetriza- 
tion effects. In the limit v — * oo, i.e. for vanishing 
phase-space density of the source, the residual correla- 
tions disappear, and so does the minimum in the cor- 
relation function. The minimum also disappears in the 
opposite limit, v — > 1, but for a different reason: in this 
limit the HBT-radii go to zero, and the whole correlation 
function becomes entirely flat. 

For multiplicity-selected event samples (as they are 
routinely collected in heavy ion collisions because 
multiplicity serves as a measure for collision centrality) 
the absence of a minimum in the correlation function, 
combined with a clear positive HBT signal at q = 0, 
can thus be taken as evidence against strong multi-boson 
symmetrization effects and provides a lower limit on the 
phase-space volume (upper limit on the phase-space den- 
sity) of the source. The significance of this limit is con- 
trolled by the statistical error on the correlation function. 
The conspicuous absence of a minimum in the correlation 
data known to us supports the evidence for low freeze-out 
phase-space densities (and thus weak multi-boson effects) 
which was obtained with other methods fl6[ ]. 

On the other hand, searching for such a minimum 
in event ensembles with approximately fixed multiplici- 
ty appears to be the most direct approach to find sig- 
natures of multi-boson symmetrization effects in heavy 
ion collisions, since it does not rely on model fits to the 
correlation function. It must be noted, however, that 
multiplicity-averaging tends to smear out the minimum, 
making it dis appear, for instance, in the "Poisson limit" 
of Sec. Ill C . To obtain a clear signature thus may re- 
quire a sharp selection on multiplicity. In practice the 



method will be limited by event-by-event fluctuations in 
the slope of the momentum spectra (A) and of the freeze- 
out radius R which both enter the phase-space volume v. 
Detailed simulations of such smearing effects go, how- 
ever, beyond the scope of the present paper and will be 
discussed elsewhere. 

Such a search for multi-boson effects via residual cor- 
relations requires a 2-Gaussian fit to the correlation func- 
tion like the one obtain ed by com binin g according to 
Eq. ( |5.]|) the expressions ( |5.12 ) and Q5.21 ). Compared to 
the usually employed fit functions this doubles the num- 
ber of HBT radius parameters and thus puts very high 
demands on event statistics. By triggering on central 
collision events in order to exploit azimuthal symmetry 
of the emission region, the number of radius parameters 
in the single-Gaussian fit can in general be reduced to 
four p2| (instead of the only two parameters required 
for the spherically symmetric toy sources studied here). 
For equal-mass collision systems, these can be further 
reduced to two parameters (R± and Rl, perpendicular 
and parallel to the beam direction) by selecting pairs with 
vanishing pair momentum in the center-of-mass system 
p2fl . To describe the residual correlations then requires 
only two additional parameters, rj_ and ri,. The com- 
plexity of the resulting fit problem is then similar to 
the now routinely applied 3-dimensional HBT analyses 
in relativistic heavy-ion physics p^| , which renders the 
proposed search practicable. 



VII. CONCLUSIONS 

We have studied multi-boson effects on two-particle 
correlations for systems with a variety of multiplicity dis- 
tributions. In contrast to the previously studied case 
of a fixed Poissonian input multiplicity distribution at 
v — ► oo, which is then at finite v modified by multi-boson 
effects, we find that in the general case the two-particle 
correlation function is affected not only by Bose-Einstein 
correlations, which can be used to extract source size in- 
formation, but also by residual correlations which depend 
on the pion multiplicity distribution, the average phase- 
space density and the pair momentum. This is true in 
particular for events with fixed multiplicity. 

These residual correlations affect the asymptotic nor- 
malization of the correlation function at infinite relative 
momentum and the "correlation strength" , i.e. the ratio 
of the correlation function at zero and infinite relative 
momentum minus 1. The asymptotic normalization de- 
pends on the multiplicity distribution, in particular the 
mean multiplicity (n), and the phase-space volume v . 
The correlation strength depends additionally on the pair 
momentum K. Residual correlation effects can lead to 
correlation strengths which are either much larger than 
1 (for example for the Bose-Einstein multiplicity distri- 
bution at intermediate values of v) or much smaller than 
one and even negative (for example for fixed multipli- 
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city in the quantum saturation limit v — > 1). All this 
happens for completely chaotic sources as those studied 
in the present paper and has nothing to do with par- 
tial phase coherence in the particle emission process or 
other quantum optical phenomena. The residual corre- 
lations arise, even in the absence of HBT correlations, 
from the generically different weights of pairs, triplets, 
quadruplets etc. of particles within a single event and 
from mixed events. The former enter in the numerator, 
the latter in the denominator of the correlation function. 
These different weights play a role when multi-particle 
symmctrization effects begin to modify the phase-space 
distribution of such pairs, triplets, quadruplets ets. of 
identical particles, i.e. for sufficiently large phase-space 
densities. The residual correlations vanish only acciden- 
tally for very specific multiplici ty dis tributions, like those 
of the "Poisson limit" in Sec. Ill C or of the "covariant 
current formalism" in Appendix [A[ 

The separation of residual and genuine Bose-Einstein 
correlations requires a fit of the correlation function with 
two 3-dimensional Gaussians. It thus complicates consid- 
erably the usual HBT analysis for extracting source sizes 
from two-particle correlation functions. Fortunately, we 
could show that the range of the residual correlations 
and of the Bose-Einstein correlations are defined by well- 
separated length scales such that, at least with suffi- 
ciently accurate data, the separation is possible. In prac- 
tice, the demands on the statistical accuracy of the data 
are formidable, and in Section VI we therefore suggested 
to begin such studies by restricting the fits to certain 
kinematic regions (small momenta in the center-of-mass 
system in collisions between equal-mass nuclei) where 
only a reduced number of fit parameters is required. 

We have also made a comprehensive model study of 
multi-boson effects on the HBT radii extracted from the 
Bose-Einstein part of the correlation function. In the ab- 
sence of residual correlations, we showed that it is possi- 
ble to isolate the source size information by first dividing 
out a contamination from ratios of single-particle spec- 
tra, using measured spectra. In the presence of resid- 
ual interactions this is no longer possible, and the HBT 
radii cannot be corr ected for this unwanted contamina- 
tion. In Section VD| we showed how to measure it within 
the model and saw that it becomes dangerously large for 
sources with high phase-space densities, i.e. when multi- 
boson effects become essential. 

It is quite likely that the present study eventually turns 
out to be rather academic since Nature only allows par- 
ticles to decouple from the collision region after it has 
had time to expand to such low phase-space densities 
that multi-boson symmetrization effects become negligi- 
ble. Existing data analyses seem to point in that direc- 
tion [jl5 16 1 . However, until we know this for sure we 
should be open to the possibility of strong multi-boson 
symmetrization phenomena and their complicating effect 
on HBT interferometry. In Section VI we have suggested 
an essentially model-independent way to search for such 
effects which requires little or no prior knowledge about 



the source size (except for that needed to perform final 
state interaction corrections): it is based on the exis- 
tence of the minimum in the (FSI corrected) two-particle 
correlation function at non-zero relative momentum q in 
events with (approximately) fixed multiplicity. If no such 
minima are found, the results from the present paper 
can be safely shelved, and one may return with strength- 
ened confidence to the conventional HBT formalism and 
source size extraction procedure based on Eq. ( 3.20] ) . 
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APPENDIX A: RELATION TO THE 
COVARIANT CURRENT FORMALISM 



The general ansatz (2.1), (2.7) for the source density 
matrix on which the present paper is b ase d is formulated 
in terms of the unnormalized states (2jq). In this Sec- 



tion we show that the density matrix of the "covariant 
current ensemble" pO| , pi[ which instead uses normalized 
states can be brought into a similar form, with one im- 
portant difference: in the covariant current ensemble and 
in the formalism used in the present paper the normaliza- 
tion of the density operator and the integration over the 
phase-space positions at which the elementary currents 
are centered are performed in opposite order. We discuss 
some relevant implications of this difference. 

The covariant current formalism is based on the den- 
sity operator [|20|,pif 



N=i \i=r J J 

where | J) is a normalized coherent state, 
| J) = \J[N;{C t , &}]> 



\J)(J\, (Al) 



-fi/2 



exp 



dpj{p) a+ 



|0) 



(A2) 



with the current J(p) from Eq. ( 2.4) and normalization 
exp(— n/2) with n given in (2.10). It corresponds to a 
specific pion multiplicity distribution p n which can be 
calculated by evaluating 



Pn 



(n\p\n) 



(A3) 
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with the states (2.8). Since 



W) = e -n/2» 



(A4) 



each coherent state itself has a Poissonian multiplicity 
distribution: 



\{n\J)f 



Thus we have 



(A5) 



(A6) 



where n depends on all integration variables. By sum- 
ming over n and using ^ n e ~"(n "/n!) = 1 and the nor- 
malization of Pn as well as (2J2), one easily checks that 
p n is normalized: 



E&> 

71=0 



(A7) 



We can decompose the coh erent state projection oper- 
ator in (Al) into the set ( |2.8| ) of unnormalized rt-particle 
states \n): 



\j){j\= E 

n.rn— 



E 



\n)(n\J)(J\m) (m\ 
(n\n)(m\m) 

n) (m| . 



V nl ml 



(A8) 



For the calculation of the n-particle spectra, which al- 
ways involve equal numbers of creation and annihilation 
operators, only the diagonal terms n = m contribute. We 
can thus replace the density operator p of the covariant 
current ensemble ( [Al] ) by 

P' = E P » (ft / d Ci p(Ci) / ^|n)<n| . (A9) 



Writing this in the form (2.1), namely 



P' = E P~ n P n ' 

n=0 



(A10) 



one finds for the density operator of the subspace corre- 
sponding to events with fixed particle multiplicity n the 
expression 



When comparing this with the ansatz (2.7) a charac- 
teristic difference becomes apparent: the covariant cur- 
rent formalism involves an additional weight function 
exp(— n[N; <f>i}]) under the integral over the positions 
and phases of the elementary source currents jo- This 
weight results from the use of normalized coherent states. 

An arbitrary weight function would at this point stop 
further analytic progress. For the specific weight function 
associated with the use of normalized coherent states, 
however, all n-particle spectra can still be calculated ana- 
lytically J20]j2l]] , in fact with much less effort than for the 
ansatz fl2.7| ). This is due to the simplifications resulting 
from the fact that coherent states are eigenstates of the 
annihilation operator. For example, the single particle 
spectrum is given by 



iVi(p) = E p tr (p' a+ap) 



(A12) 



= E P jtPN (flj d ^(Ci) / ^) \J(P)\ 2 (J\J) 



N=l \i=l 



We write 



(J|J) 



E 



(A13) 



and calculate the normalization of the single particle 
spectrum as 



dpN^p) 



N 



EMIT / ^p(o) 



N=l \i=l 
00 / N 



En" 
— - e 
_,, , nl 

/ n=0 



N=l \i=l J J / m=0 



t in 



n" 
m — - e 
ml 



E m = ( m ) 

m=0 



(A14) 



The single particle spectrum is thus correctly normalized 
to the average particle multiplicity, calculated from the 
multiplicity distribution p n of the covariant current en- 
semble. Similarly one finds that the two-particle spec- 
trum 

N 2 ( Pl ,p 2 ) = E P1 E P2 E Pn [ II / d <iP«i) I 

x\J( Pl )\ 2 \J(p 2 )\ 2 (J\J) (A15) 



is correctly normalized to twice the average number of 
pairs calculated from p n : 



5 » = -F- E P " (li I KC<) / ^] e~ n \n)(n\ . , » 

n] P n N=i \ t =i J J 2n J dp 1 dp2N2( Pl ,p 2 ) = J2 m ( m - 1 )Pm = ( , m(m-l)). 

(All) ' m=0 

(A16) 
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From t he stu dies in (2^j2l[] it is known that the normali- 
zation ( A14 ) can be alternatively written as 

(m) = ((AT)) n , (A17) 

where ((N)} is the average nu mber of elementary currents 
and no is their normalization (2.6). The two-particle cor- 
relation function is in the covariant current formalism 
given by 

«JV(JV-1)» (, , \J x S(x,K)e^\ 2 \ 



C(q,K) 



amy 



i 



J x S{x, Pl )JS{x,p 2 ) 



(A18) 

where the effective emission function is nothing but the 
Wigncr density associated with the currents J: 

d *y c -iKy 



S{x,K) 



2(2tt) 3 



(A19) 

This is, up to the prefactor, the standard form ( p. 20 ). 
The covariant current ensemble with its multiplicity dis- 
tribution p n is thus a secon d example (besides the "Pois- 
son limit" of Section III C) for a system in which multi- 
boson effects do not cause residual correlations. Up to 
the different normalization (((N(N— 1)) ) vs . 1), the only 
difference between the density operator (2/7) in the "Pois- 
son limit" and th e c ovariant current ensemble with the 
density operator ( |A9| ) is the different definition of the ef- 
fective emission function: in the former case it is given 



by Eqs. ( 3.18 ) and ( 3.19 ) while for the latter the much 
simpler definition (A19), with all multi-boson effects al- 
ready fully accounted for. The price to be paid is a more 
c omp licate d ex pression for the multiplicity distribution 
((|A6|) vs. ( [3.11 ) - while, for Gaussian sources, the latter 
is known analytically, the former seems to require numer- 
ical evaluation). In both cases the measured multiplicity 
distribution (p n andp n , respectively) depends in a com- 
plicated way on the mean phase-space density d of the 
source. 



APPENDIX B: CALCULATING THE 
CORRELATION FUNCTIONS 



as well as 



cosh 





f 1 


~)1 


V 




a jJ . 



(B2) 



A few simple steps, using eqs. (4.4) and (5.2) and exploit- 
ing the symmetry of the sums under exchange of i and j, 
lead from (5.1) to 



C. 



R 2 (q,n) 



YT 1= iA lJ {K)E lJ {q,K) 



(B3) 
(B4) 



For very d ilut e (v — > o o) and very dense (v — ► 1) 
sources Eqs. (4.7) and ( |4.8| ) can be used to calculate ana- 
lytically the behaviour of the correlation function in the 
limits q — > and q — > oo. Note that in these limits the 
depende nce on the angle between q and k (which en- 
ters via (B2)) disappears such that all following results 
depend only on the modulus k of the pair momentum: 

1. q — » oo : 

For v — > oo the g 2 -dependent exponential terms in the 
numerator of i?2 are much smaller than those in the de- 
nominator, and the Bose-Einstein correlations thus va- 
nish as q — ► oo. For v — > 1 the analysis is only slightly 
harder: writing v = 1 + e, one finds that the exponen- 
tial terms of both the numerator and denominator are 
given by exp[— q 2 (l — e)] for all values of i,j except for 
the lowest order term with i = j = 1; for the latter the 
numerator goes as exp[— q 2 ] while the denominator goes 
as exp[— <7 2 (1 — 2e)]. Again the numerator vanishes more 
quickly, and the Bose-Einstein correlations disappear at 
q — > oo, as they should. 



For the residual correlations (B3) one checks similarly 
that in both limits (v — ► oo and v — ► 1) for q — > oo 
numerator and denominator arc both dominated by the 
term i = j = 1, giving 



lim C\ >{q,K) = — = -2 =N 



ho 



Bn 



h\ 



(B5) 



In this Section we give some technical steps for the 
analytical evaluation of the correlation functions for the 



Gaussian sources defined in Sec. [V 
Let us define 



A-ij (^) — ^i+i c i Cj & 

-E^ij (^) — — h% hj Ci Cj c 



A{k) 
B(k) 



oo 
i J = l 



(Bla) 
(Bib) 
(Blc) 

(Bid) 



2. g 



: 



I n t his limit R2 approaches trivially the value 1 (see 
(|5.l|)), while the residual correlations are given by 



lim C 



B{n) 



(B6) 



The incoherence parameter defined in Eq. (3.22) is th us 
Aincoh(«) = 1 , whereas the correlation strength ( 3.21 ) is 
given by (|]|). 

Other useful limiting expressions are 
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lim hi = (n) , 

v — >oo 

lim h,2 — (n(n—i)) , 

v— too 

lim hi = (n(n-l) ■ • ■ (n-i+1)) for i > 2 , (B7) 

V— >oo 

lim A(/t) = lim A u (k) = (n(n-l)) e^ 2 , (B8) 
lim = lim B u (k) = {n) 2 e - K2 , (B9) 

V— >00 V — 'OO 

and (to linear accuracy in e) 

oo 

lim hi = (l + fie) ^ p n - %m for i > 1 , (BIO) 

(Bll) 



lim = M"-i)) e -K» 



1 C 2 , (n) - (l-j?o) 2 „\ ' 

1 — (k — 3)e H — — - — (2k — 6)e 



(n(n-l)) 



lim -B(k) = (n) e 



(n) - (1-po) . 2 ' 
^— r (k — 3)e 



(B12) 



These last identities require the evaluation of the follow- 
ing sums involving hi = Y^ =l Pn- 



= (n) - (1-po) , 
(n(n-l)) 



i=2 
oc 



^i+j - — 2 — ' 

~ - _ ((n-l)(n-2)) 

2 , "i+i+i 2 Po ' 



(B13) 
(B14) 
(B15) 



which is achieved by appropriate reordering and rela- 
belling of the nested sums. 



APPENDIX C: CALCULATING THE RADIUS 
PARAMETERS 



Inserting ( B4 ) into Q5.14 ) and separating the constant 
and ^-dependent terms one finds 



R 2 



id 

?2^.\ _ E>2 



Aij (k) 





1 1 




V 2 


KOLi 


Oj/J 



-.Ri(«;) 

E 2 



9 OO 

ft \ - 



1 



i; 2 ^ 2A( 



k) V o,; a,- 



(CI) 
(C2) 



In the dilute gas limit v — * oo only the lowest order terms 
i=j=l cont ribute to the sums, and one finds the expres- 
sions ( 5.15 ). In the opposite limit v — > 1 one can use 
([BTol)-(|B12D to obtain 



lim fi 



(1-Po) 



(n(n — 1)) 

2 



(C3a) 



, (ra)-l+e~W 



for fixed n > 2 , 
Bose— Einstein , 

Poisson , 



lim r2 = . 

U->1 



(C3b) 



All these limits are K-independent; however, in the in- 
termediate region 1 < v < oo the HBT radii develop 
a K-dependence as shown by the numerical results in 
Sec.[VD. 

The vanishing of r 2 in the limit v — > 1 has an inter- 
esting origin: in this limit the functions Gi(p 1 ,p 2 ) in 
(4.4) beco me all identical, and the single particle spec- 
trum ( |3.8| ) is again a simple Gaussian with width A, just 
as in the opposite limit v — > oo. In this case the angle- 
dependent term ~ q-K drops out from the product of 
single particle spectra in the denominator of the Bose- 
Einstein correlator R 2 , and the two HBT radii R± and 
Rh become equal. 

The "residual corre lation radii" are obtained by in- 
serting (B3) into ( 5.22 ) and separating the constant and 
^-dependent pieces: 



A(k) 



R 2 NB{k)-A{k) 

'Bij(K) Aij(K) 



(C4) 



1 \ - ( Bijjn 
■ 2v 2 .4^ V B(k) 



?(k) A(k) 



R 2 

,2 



2u 2 ( 



i,3'=l 



Af B(k) - A(k) 
Aj(k) _ Bjj (k) 
A(k) B(k) 



1 N 2 
ft, 



(C5) 



To facilitate the discussion in the main text (Section |VE| ) 
we here give a specific e xpre ssion in the limit v — > 1. We 
decompose the r.h.s. of (C4) into three factors which we 
expand in powers of e = v — 1: 



lim ^4(k 

ti-»l+e 



a + a! e + . . . 

lim (NB(k) -A(k)) =b 

1) 2 ' 



he- 



lim 

u->l+e 2w 



c + ci e - 



(C6a) 
(C6b) 

(C6c) 



The lowest order terms have been calculated (see also 

mmy. 



«0 

&0 
CO 
Cl 



(n(n-l)) 



1-PO-Pl/ \2 



i(n-l)) 



(1-Po) 



0, 



1 ~Po o (n(n-l)) 
' (n) 2(n} 2 



(C7a) 

(C7b) 
(C7c) 
(C7d) 
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Using Eq. (C3a) in the form 
lim R 2 ± (k) = 2 



one finds 



-' o( n ) _ O-PO) 2eR 2 



n(n— 1)) 



,. r±(K) _ 1 (n(n-l)) a ci 
£Sitf («) -4(n)-(l-p ) 6 



(C8) 



(C9) 



This relation is the basis of the discussion in Sec. VE 



For the Bose- Einstein multiplicity distribution (5.3d) 
an unfortunate coincidence causes bo th bo and c\ to va- 
nish. In this case their ratio in ( p9| ) must be replaced 
by Oz/b\. The calculation of C2 requires a consistent cal- 
culation to order e 2 which we have not done. Therefore, 
for the case of a Bose-Einstein distribution we have no 



analytical expression for the limit (C9) 



[1] W.A. Zajc, Phys. Rev. D 35, 3396 (1987). 

[2] S. Pratt, Phys. Lett. B 301, 159 (1993); Phys. Rev. C 

50, 469 (1994); S. Pratt and V. Zelevinsky, Phys. Rev. 

Lett. 72, 816 (1994). 
[3] W.Q. Chao, C.S. Gao, and Q.H. Zhang, J. Phys. G 21, 

847 (1995); Q.H. Zhang, W.Q. Chao, and C.S. Gao, Phys. 

Rev. C 52, 2064 (1995). 
[4] A. Bialas and A. Krzywicki, Phys. Lett. B 354, 134 

(1995); A. Bialas and K. Zalewski, Phys. Lett. B 436, 

153 (1998); and Eur. Phys. J. C 6, 349 (1999). 
[5] J. Wosiek, Phys. Lett. B 399, 130 (1997); K. Fialkowski, 

R. Wit and J. Wosiek, Phys. Rev. D 58, 094013 (1998). 
[6] Q.H. Zhang, Phys. Lett. B 406, 366 (1997); Phys. Rev. C 

57, 877 (1998); Phys. Rev. C 58, R18 (1998); and Nucl. 

Phys. A634, 190 (1998). 
[7] T. Csorgo and J. Zimanyi, Phys. Rev. Lett. 80, 916 

(1998) ; J. Zimanyi and T. Csorgo, Heavy Ion Physics 
9, 241 (1999). 

[8] UA. Wiedemann, Phys. Rev. C 57, 3324 (1998). 
[9] R.L. Ray, Phys. Rev. C 57, 2523 (1998); R.L. Ray and 
G.W. Hoffmann, Phys. Rev. C 60, 014906 (1999). 
[10] Q.H. Zhang, P. Scotto, and U. Heinz, Phys. Rev. C 58, 
3757 (1998). 

[11] S.V. Akkelin and Yu.M. Sinyukov, Nucl. Phys. A661, 
613c (1999). 

[12] A. Bialas and K. Zalewski, Phys. Rev. C 59, 097502 

(1999) . 

[13] R. Lednicky et al., Phys. Rev. C 61, 034901 (2000). 
[14] G.F. Bertsch, Phys. Rev. Lett. 72, 2349 (1994). 
[15] J. Barrette et al. (E877 Coll.), Phys. Rev. Lett. 78, 2916 
(1997). 

[16] D. Ferenc, U. Heinz, B. Tomasik, U.A. Wiedemann, and 
J.G. Cramer, Phys. Lett. B 457, 347 (1999); B. Tom asik, 
U. Heinz, and U.A. Wiedem ann, tiucl-th/9907096 . 

[IT 



J. Rafelski and J. Letessier, iucl-th/9903018 



[18] Q.H. Zhang, Phys. Rev. C 59, 1646 (1999). 

[19] We will repeatedly use the phrase "Poissonian input mul- 
tiplicity distribution". In order to avoid confusion we 
stress that by this we mean that the multiplicity dis- 
tribution is Poissonian in the infinite phase-space vol- 
ume limit v — ♦ oo, but differs from a Poissonian at fi- 
nite values of v in a specific, n-dependent way (see for 
example Eq. (3.11)). It does not mean that at finite v 



J. Letes- 



sier, A. Tounsi, and J. Rafelski, nucl-th/991104; ; J. Le- 
tessier and J. Rafelski, micl-th/0003014[ 



[20 
[21 

[22 
[23 



[24 
[25 

[26 



[27 

[28 
[29 



[30 



[31 

[32 
[33 



[34 

[35 

[36 
[37 



one starts out with a Poissonian multiplicity distribution 
and that multiboson effects then change this input into 
something else. Multiboson symmetrization itself cannot 
change the multiplicity distribution, as imprecise formu- 
lations (e.g. in Q) may have sometimes suggested: start- 
ing from exactly n classical (distinguishable) particles in 
a given phase-space volume and then symmetrizing the 
multiparticle wave function still leaves us with exactly n 
particles, not more and not less. 

M. Gyulassy, S.K. Kauffmann, and L.W. Wilson, Phys. 
Rev. C 20, 2267 (1979). 

S. Chapman and U. Heinz, Phys. Lett. B 340, 250 (1994); 
U. Heinz, in Correlations and Clustering Phenomena in 
Subatomic Physic (eds. M.N. Harakeh, J.H. Koch, and O. 
Scholten), N ATO ASI Series B 359, 137 (1997) (Plenum, 
New York) fcucl-th/9609029| ). 

U.A. Wiedemann and U. Heinz, Phys. Rep. 319, 145 
(1999); U. Heinz an d B.V. Jacak, Ann . Rev. Nucl. Part. 
Sci. 49, 529 (1999) flnucl-th/990202C| ). 
W.A. Zajc, in Particle Production in Highly Excited Mat- 
ter (eds. H.H. Gutbrod and J. Rafelski), NATO ASI Se- 
ries B 303, 435 (1993) (Plenum, New York). 
C. Slotta and U. Heinz, Phys. Lett. B 391, 469 (1997). 
P. Scotto, PhD thesis, University of Regensburg, Sept. 
1999; P. Scotto and U. Heinz, manuscript in preparation. 
W.Q. Chao, C.S. Gao, and Q.H. Zhang, Phys. Rev. C 49, 
3324 (1994). (In Eq. (2) of that paper the normalization 
factor was inadvertently omitted but the rest of the paper 
shows that normalized coherent states were used.) 
C. Slotta, PhD thesis, University of Regensburg, Jan. 
1999. 

E.V. Shuryak, Phys. Lett. B 44, 387 (1973). 

R.M. Weiner, Phys. Lett. B 232, 278 (1989); 242, 547 

(1990); M. Plumer, L.V. Razumov, and R.M. Weiner, 

ibid. 286, 335 (1992); I.V. Andreev, M. Plumer, and 

R.M. Weiner, Int. J. Mod. Phys. A 8, 4577 (1993). 

M. Biyajima, A. Bartl, T. Mizoguchi, N. Suzuki, and O. 

Terazawa, Prog. Theor. Phys. 84, 931 (1990); N. Suzuki 

and M. Biyajima, ibid. 88, 609 (1992). 

J.G. Cramer, Phys. Rev. C 43, 2798 (1991); J.G. Cramer 

and K. Kadija, ibid. 53, 908 (1996). 

U. Heinz and Q.H. Zhang, Phys. Rev. C 56, 426 (1997). 
U. Heinz, in Measuring the Size of Things in the Uni- 
verse: HBT Inter ferometry and Heavy Ion Physics (eds. 
S. Costa, S. Albergo, A. Insolia, a nd C. Tuve), p. 66 
(World Scientific, Singapore, 1999) ( |hep-ph/9806512[ ). 
U.A. Wiedemann, D. Ferenc, and U. Heinz, Phys. Lett. 
B 449, 347 (1999). 

H. Merlitz and D. Pelte, Phys. Lett. B 415, 411 (1997); 
Z. Phys. A 357, 175 (1997). 

U.A. Wiedemann et al, Phys. Rev. C 56, R614 (1997). 
K. Geiger, J. Ellis, U. Heinz, and U.A. Wiedemann, Phys. 



21 



[38] 



Rev. D 61, 054002 (2000). 

The astute reader will have noticed that the spec tra and 



[39] 
[40] 



[41] 
[42] 

[43] 



QIC 



do de- 



correlations in the "Poisson limit" of Sec. 
pend on no. The reason is that in this case the measured 



multipl icity distribution p n given in (3.11) depends on 
no via (3.12). However, it really depends only on the pro- 
duct noN — by rescaling the normalization no of the 



elementary currents and the parameter TV in (3.11) si- 



multaneously, keeping their product fixed, one can either 

set no = 1 or N = 1. In other words, the normalization 

no of the elementary currents does not enter any physical 

results directly, but only via the multiplicity distribution 

if the latter is allowed to explicitly depend on no- 

U.A. Wiedemann , J. Ellis, U. Hein z, and K. Geiger, in 

RefJpJ, p. 135 ( |nucl-th/9808043D . 

In [pTfTEq. (38), an incorrect expression was given for R% 

which resulted from an incorrect rewriting of the correct 

expression given in |^|. Since the numerical simulations 

in JS^] were compared with the original correct analytical 

expression, all figures in j37| are correct, although the 

discussion in the paragraph following Eq. (38) in [^ is 

inaccurate and slightly misleading. 

W.A. Zajc et ai, Phys. Rev. C 29, 2173 (1984). 

S. Chapman, P. Scotto, and U. Heinz, Heavy Ion Physics 

1, 1 (1995). 

S. Pratt, Phys. Rev. C 56, 1095 (1997). 



22 



